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Self-gravitating systems have been investigated in various situations. Boson stars ( [[J 
for reviews) have a very simple constituent, a complex scalar field which is bound by gravi- 
tational attraction. Thus the boson star provides us with the simplest model of relativistic 
stars. 

The solutions for relativistic boson stars are only numerically obtained in four dimensions. 
In (2 + 1) dimensions, static equilibrium configurations have been argued fU in Einstein 
gravity with a negative cosmological constant. 

In this paper, we obtain an exact solution for a nonrotating boson star in (2 + 1) dimen- 
sional gravity with a negative cosmological constant. We consider that the scalar field has a 
strong self-interaction. An infinitely large self-interaction term in the model leads to much 
simplications as in the (3 + 1) dimensional case ||. The study of exact solutions will lead to 
a new aspect of gravitating systems and clarify the similarity and/or the difference among 
the other dimensional cases. 

We consider a complex scalar field with mass m and a quartic self-coupling constant A. 
The action for the scalar field coupled to gravity can be written down as 



S = Sx\f^g 



1 ,„ x ,_ ,2 2 I |2 ^ 



(R+2C) - |V^|-m 2 M' 



(1) 
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where R is the scalar curvature and the positive constant C stands for the (negative) cos- 
mological constant. G is the Newton constant. 

Varying the action ([!]) with respect to the scalar field and the metric yields equations of 
motion. The equation of motion for the scalar field is 

VV~ m 2 (p- \\(p\ 2 (p = 0, (2) 

while the Einstein equation is 
Rnv ~ ^9^R = 8vrG 2Re {v^Vyip - ~ | V<^| 2 - m 2 |<^| 2 - ^ |y3| 4 g^ + Cg^. 



(3) 

We assume the three-dimensional metric for a static circularly symmetric spacetime as 



ds 2 = -a 2 (r)dt 2 + f3 2 (r)dr 2 + ^V-dO 2 
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(4) 



where a, /3, and 7 are functions of the radial coordinate r only, while 70 is a constant. 
Because of the coordinate invariance with respect to r, we have a freedom in choosing the 
function (3(r). We will use this residual "gauge choice" later. 

We also assume that the complex scalar has a phase which is linear in the temporal 
coordinate: 



v = e -^V(r) 



(5) 



where <p(r) is a function of the radial coordinate r only and u is a constant. 
Then the equations @ and (H) are written as 

1 d I cry dip^ 



u 
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a(3 dr \f3dr) 
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1 { d<p\ uj 2 2 



A 



+ C, 
+ C, 

+ c. 



For convenience, we rescale the variables as 

A 



r = mr, <p 



8nGm 2 



tu = uj/m, and C = C/m 2 



(6) 
(7) 

(8) 

(9) 



(10) 



Using these variables, we can rewrite the equations as 



u 2 



1 d ^^_~_ A ~3 = 0; 



a 2 a/^7 dr \ P dr 



1 d (ld^ 

P'y dr 1/3 c?f 



~x 2 



~,2 



A 



/? 2 I df i a 2 ^ * 2^ ' 



1 c?a o?7 1 I dip 



_\ 2 



~,2 



+ - 



A 



a/3 2 7 <ir <ir /3 2 \ dr 



a" 



1 c/ /l da\ 1 (dip\ 2 Cj 2 _2 ~ 2 ^ _ 4 



(11) 
(12) 

(13) 

(14) 
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Next, we rescale the variables again, to study the limit of large self-inteaction. New 
variables are: 



/VA, (p* = VXip, and = AC. (15) 



= r 

Then the equations can be written as 

1«' 7' 11,2 ^ 2 9 2 1 4 

1 /a'\ 1 1 . ,. 2 £> 2 9 9 1 A _ . , 

where ' denotes the derivative with respect to r*. 

For the limit of large self-coupling, A — > 00, these equations will be reduced to[| 

Co 2 

<p*-<p*-<pl = 0, (20) 



a 2 

/? 2 a 7 
1 /W 



V* - = 


= 0, 




UJ 2 2 
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Now, let us solve the set of an algebraic equation and differential equations ( p0] - |2~3|) . 
From Eqs. (0) and fl2|), we find 



(21) 
(22) 
(23) 



a' 



0. (24) 



To solve this, we use the residual gauge choice on (3. For simplicity, we take 



1 For a finite C*, the actual value of C = m 2 C*/A becomes infinitely small if the limit A = 00 
is taken literally. We can however interpret the limit as an approximation of a large self-coupling, 
and A is not simply taken as a mathematical infinity. 



4 



a/?7 = n- (25) 

Then Eq. fl2H]) is solved as 

a 2 (r)=A(r 2 + r 2 ), (26) 

where A and r are integration constants. Since the constant scale factor of a can be 

absorbed in the rescaling of t, we can fix the value A = 1. Then, for the later use, we notice 
that 

S- L (27) 

From Eqs. (0), (0) and @, we fincf] 



77' 1 i 



2 ^ + C*. (28) 



Here using Eq. fl2~T|), we obtain 



c/ 7 2 1 



da 2 2 

We have to solve Eq. fl2"U| ) by using the solution for Eq. fl2"0| ) can be solved easily. 
The interior solution for ip* is 

& = ~ 2 ~ 1- (30) 

or 

This describes the configuration of the boson field inside a boson star. On the other hand, 
outside the boson star, 

= 0. (31) 

Before solving the interior and exterior solution for the metric, we rewrite the radial line 
element as 



2 This equation (|28|) is consistent with the other equations (20) and (21) 
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V mr V m ocy "7 



A 1 



m 2 a i^t + C* 



(32) 



Consequently, we get 



I3 2 dr 2 



A 1 



m 2 a 2 



T dY 



(33) 



First, we solve the exterior solution. Since = outside the boson star, the solution 
of Eq. (E3) is found to be 



o? = irl 2 ~ B, 



(34) 



where B is an integration constant. 

The full line element is obtained, by substituting Eqs. (p5f), (|33| ) and (j34|) into (f|), as 

1 



2 _ 1 />r »,2 ^2 D \ 7.2 ^ 



rfs -~ni\p*T-CiB dt z + , 

C 2 v ' m 2 C7 2 - C 2 5 75 



d 7 2 + \d6 2 . 



(35) 



Here we must remember that C* = CA/m 2 . We further introduce a new radial coordinate 
R defined as 



R 




7, 



(36) 



and set 7 = Jm 2 /A. Then we find 



di? 2 + flW. 



(37) 



After rescaling ^k-dt — > c?i, we find the metric is precisely the same as the well-known BTZ 
vacuum solution || . Thus we identify the BTZ mass M of the boson star with 



8GM n = CfB. 



(38) 



Here we turn to solving the interior solution. From Eqs. (|29"D and (|30|), we obtain a 
differential equation: 

2 U 2 ) 



d 7 2 



do? 



C* ■ 



(39) 
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The solution of this equation is given by 

7 2 = -~ - ^ In t 2 + \o? + C. (« 2 + D) , (40) 

where D is a constant. The value of D is determined from the condition at the boundary 
of the boson star, where if* = 0. One can find the condition for a smooth connection of the 
exterior and interior solutions (from Eqs. (p3|), ( [54]) and (|38|)):P| 



(41) 



From Eqs. (]4lD and (||]), we get 

C 2 D = 8GM . (42) 



What does the metric at the center of the boson star look like? As seen in Eq. (p7|). 
the origin is located at R = 0, or 7 = 0. Therefore «o = «(7 = 0) satisfies the following 
equation: 

- \- 2 - Cu 2 In % + \a\ + C„ag + -^8GM G = 0. (43) 
Thus the metric at the center of the boson star is written by using a as 



1 1 



ds 2 (R ss 0) sa -a^t 2 + —7— ^rf/? 2 + (44) 



where (y9 2 = ^ — 1 is the value of at the center of the star. We can give an arbitrary 
value for (p 2 by choosing the value of uj. Thus we can take ip 2 Q as a free parameter which 
characterizes a solution for a boson star, instead of Co. 

In higher dimension more than three, we can fix the natural value of the metric compo- 
nent at the center. In three dimensions, we cannot take a definite choice, because there is 
no local definition of mass; an additional constant is allowed because it does not lead to a 
curvature singularity directly, unlike the higher-dimensional cases. This is also the reason 



3 The first derivative is automatically connected when this condition (|4l| ) is satisfied. 
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why an arbitrary constant remains in the study of Cruz and Zanelli on relativistic stars 
in (2 + 1) dimensions. 

Nevertheless, we can choose a plausible choice of the metric at the center,^ in order to 
give the definite properties of the boson star; if one wishes to know the general case, one 
need only loose the constraint in the following analysis. 

We choose the following condition: 

agQvHb + C.) 2 = l- (45) 

According to this choice, the spatial metric at the origin looks like the one with no conical 
singularity. The condition fl4"5|) determines the value of « , as a function of v?*o- 

This choice for the boundary condition makes the mass of the scalar boson definite. In 
other words, if we take other values, the parameter m does not represent the mass of the 
excitation of the boson field p near the origin; the value must be rescaled. In the present case, 
since there is no asymptotically flat spacetime region, the identification of mass parameter 
in the most interior region is reasonable. 

We obtain the relation between the BTZ mass of the boson star and ip*o from Eqs. (f43|) 
and g|): 

8GM Q 1 



\<P% + flo ~ (flo + l) In (<^o + l) - C* 



(46) 



The value of BTZ mass M a is negative for small <£>*o- in particular, as (p^Q goes to zero, 
M a approaches —1/(8G). One can take 

Mi = M Q + ^, (47) 

as a definition of the intrinsic mass M$, we call here. When ip^ vanishes, Mj becomes zero, 
and spacetime becomes a usual anti-de Sitter spacetime everywhere, i.e., 



4 The conical sigularity at the origin is not a priori governed by the differential equations because 
only the boundary condition is concerned. This is even true for a vacuum case. 
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ds 2 = - (l + CR 2 ) dt 2 + dK + R 2 d6 2 , (48) 

where the temporal coordinate has been appropriately normalized. Then we have so-called 
AdS vacuum.^ 

The total particle number N is given by 

N = J d 2 x^g \g tt \ i (<pd t tp* - <f*d t if) . (49) 

Therefore we can calculate the particle number of the boson star by making use of the above 
solution. It is given by 

AT O f 2 1 a Pl A 

J 87rGA a 2 70 

mCo r (Co 2 \ I Ha~ , 
— - 1 — n\ — dr* 



2GA70 J \a 2 J a 2 V m 2 
mCo I A r^ 2 I 'Co 2 \ 1 



1 —da 



2 



4GA70 V m 2 Joel \a 2 J a 2 

m I A Co fed 2 1 , Co 2 " 
r"o — I 9 1 hi 



4GA70 V fn 2 OiQ \ Oq «o 



1 Jcplo + 1 [^o - In (cpl + 1)] 



AGm y%+ a ' (50) 

The intrinsic mass and particle number of the boson star as functions of (p^ are plotted in 
FIG. § 

For small Lp*o, the mass and the partcle number are given approximately by 

M ^ mN ~ i&k; (51) 

Both maxima of the mass and particle number as functions of <£>*o are located at ip*o = 
f*om- f*om satisfies the following equation: 

6C^l 0m - 4^ 0m ~ V% m + (-2C, + 4^ 0m + 3^ 0m ) In (l + rf 0m ) = 0. (52) 



5 Note that, in the presence of matter, the limit of "no matter" cannot yields so-called black hole 
vacuum, which is given by BTZ geometry in the limit M Q = |]. 
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Unfortunately, the solution of this equation cannot be expressed by the usual mathematical 
functions. For small C*, we can solve the equation (|52| ) approximately and get y?*o m ~ 
(24C,) 1 /6. 

The maximum values of and N, which are denoted by Mi max and N max respectively, 
is shown in FIG. |2] as functions of C*. It is worth noting that the value of maximum BTZ 
mass (M D max = Mi max — ^) is positive for all finite values of C*. In particular, in the 
limit C* = 0, the value of M a max vanishes. The behavior of M Q max for small C* is given by 
M max « a 2/3 /(32^3G). The value of N max in the limit C* = is l/(4Gm). 

The binding energy Ej, max = Mi max — mN max for the maximum case is negative for 
sufficiently small C*, as seen from FIG. |2|. This implies that the boson star solution for 
small C* is energetically stable. When the value of C* is larger than a critical value C* cr u, 
Eb max turns to be positive. We find that the critical value C* cr i t is 

243 io + 4v /l3-31n^ 

To see the meaning of the critical value, we plot the dependence of binding energy 
Eb = Mi — mN on the particle number N for a fixed C*. FIG. [3] shows the so-called 
bifurcation diagrams f|. The cuspoidal point corresponds to the boson star with the mass 
Mi max and the particle number N max . This point represents an absolute stable solution ||. 
When C* > C* cr u, however, the cuspoidal point is not a global minimum of the binding 
energy. In this case, the boson star is not stable because the binding energy is positive for 
all finite iV.Q 

Lastly, we discuss the size of the boson star. We define the radius of the boson star by 
the radius of the boundary of the interior and exterior solutions, where y?* = 0. The radius 
of the boson star Rb is derived from Eqs. (HO) and (|4"2"D with the definition (|36|). We find 



16 1598 + 434^ - f 1242 + 351^) In ■ 
C; crit = ^_ —1 5_ « 2.5268. (53) 



6 Notice also that the point corresponding to N = and Eb = is always a minimum point, 
because Eb ~ +(^ /(24GC*) for small ^*o- 
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2 A 

R b = — 

mr 



A 



C»(^ 2 o+l) 8GM 

, 2 + ^ 



1 



771' 



T^o + flo - (^*o + l) In (<^o + l) + C*^q 



(54) 



The radius of the boson star with maximum mass as a function of C* is plotted in FIG. 



For small C*, i4 of the boson star with maximum mass behaves as Rh ~ a/ -hf — t7b ■ 
The profile of the boson field can also be solved in the same way. That is given by 



R 2 



A 



2 (^*o + 1 



2 + C* 
^+1 



V*o + 1 . 



■ (55) 



Of course R^ ([>|) is also obtained from this equation by setting ip* = 0. One can see from 
(|55|) that a has no singular behavior, because a 2 = uj 2 / (ip 2 + 1). The profiles of the boson 
field if 2 for specific values of C* are shown in FIG. |5|. 

To summarize, we have obtained an exact solution describing a boson star with very large 
self-coupling constant in (2+1) dimensions. There is a critical value for C* = (A/ (87rGm 4 ))C, 
C* crit ~ 2.5268. For C* > C* crit the binding energy turns to be positive, the boson star 
configuration is not energetically favorable. 

Explicit study of stability of the boson star under linear perturbations and pulsations 
will be discussed in future publications. The analysis of the boson star with an arbitrary, 
actually finite value of the self-coupling is of much interest and will be necessary. 

The rotation of the boson starf] is easily incorporated in the (2 + 1) dimensional model. 
The rigorous result on a nonrotating boson star in this paper is a basis of the study of the 
rotatinge case. 

More general cases including such as a \(p\ 6 coupling may exhibit more complicated 



results, but the analysis of them can be carried out in the same manner as in this paper. We 



7 Properties of a spinning boson star with large self-interaction in (3 + 1) dimensions is studied by 
Ryan §. 
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expect that our analysis can also be easily extended to the solutions for fermion-boson stars 
in three dimensions. It will be the subject to be studied elsewhere. 
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FIGURES 
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FIG. 1. Boson star mass Mj in units of G 1 (solid line) and particle number N in units of 
(Gm)" 1 (broken line) as functions of ip*o for different values of C*: (a) C* = 1. (b) C* = 2.5268. 
(c) C, = 4. 
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FIG. 2. The value of maximum intrinsic mass Mi rnax in units of G 1 (solid line) and The value 
of maximum particle number N max in units of (Gm)~ l (broken line) as functions of C*. 
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FIG. 3. Binding energy Ej, (in units of G -1 ) as a function of N (in units of (Gm)^ 1 ) for C* = 1 
(solid line), C* cr it(~ 2.5268) (broken line), and 4 (dotted line). 
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FIG. 4. The radius of the boson star i?5 (in units of y/ 'A/m' 2 ) with maximum mass Mj max as 
a function of C*. 
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FIG. 5. The behavior of the scalar field (p 2 as a function of R (in units of ^ A/m 2 ) for a boson 
star with maximum mass Mj ma:E for C* = 1 (solid line), C* C r-it(~ 2.5268) (broken line), and 4 
(dotted line). 
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